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Abstract 

Surveys of redshifted emission lines in the near-infrared will provide unique information 
on the cosmic star formation history. Near-infrared Ha surveys will probe the cosmic star 
formation density and the evolution of the luminosity function of star forming galaxies 
out to z ~ 2.5, including the important z = 1.5 — 2.0 epoch, where according to current 
knowledge the cosmic star formation density peaks. Ha and other hydrogen recombination 
lines (except Lya) and [O ii] are all useful tracers of star formation provided the luminosity 
function is probed to sufficiently deep levels. However, shallow surveys or surveys of Lya 
or [Oiii] emission are not suitable for measuring or even constraining the cosmic star 
formation density. A first estimate of the star formation density at z ~ 2.25 is derived, 
based on a new and sensitive large-area ESO near-infrared survey for Ha emission at this 
redshift. Future observational directions, at near-infrared and longer wavelengths, for 
reading the cosmic star formation history are outlined. 


1 Introduction 


The very rapid developments in near-infrared (near-IR) array technology of the last few years 
have brought about a revolution in the held of emission line surveys in the near-IR wavebands. 
Until recently such surveys did not yield any detections |^, 1^,0]! although these null-results in 
some cases began to set interesting limits on the number density of the most luminous emission 
line galaxies |^. The recent availability of near-IR cameras with large-format arrays has now 
produced the hrst detections in surveys for Ha emission at high redshift. While the number of 


Invited Review to apper in “Extragalactic Astronomy in the Infrared”, eds. G. A. Mamon, Trinh Xuan 
Thuan, & J. Tran Thanh Van, Editions FrontiHes, Gif-sur-Yvette 


1 




detections is still small, this result for the hrst time provides a crude estimate of the density 
of luminous emission line objects at cosmological distances. Simultaneously, a change of focus 
has taken place. Whereas previously surveys were aimed at finding only the most luminous 
objects, current technology allows deeper searches, and for the first time an estimate of the 
star formation density (SFD) at z ~ 2.25 is obtained. Future surveys will be able to probe the 
luminosity function to even greater depth. 

This review first discusses our current state of knowledge of the star formation history of 
the universe, and the crucial role which near-IR emission line surveys will play in this field 
(Sect. 1^). Subsequently the use (and limitations) of emission lines as cosmic star formation 
tracers is discussed (Sect. Then, results of the ESO near-IR Ha survey are presented, and 
implications for the cosmic SFD at z ~ 2.25 are derived (Sect. Finally, prospects for future 
observational work, using Ha and other tracers, for tracing the star formation history of the 
universe are discussed (Sect. H). Throughout this paper the Hubble constant is written as 
Hq = 75/175 km s“^ Mpc“^ and go = 0.1 is assumed. 


2 Cosmic star formation history 


The star formation properties of the universe at high redshift form one of the focal points of 
current observational research in cosmology. Recently, our current understanding of the cosmic 
star formation history has been summarized 0 by combining the local Ha objective prism 
survey results p and low redshift {z < 1) Canada-France redshift survey (CFRS) results |]T4 
with estimates of the SFD at 2 ; ~ 2.75 and ^ ~ 4 based on the Hubble Deep Field (HDF) data. 
These combined data sets indicate a SFD rapidly increasing with cosmic time at high redshift, 
reaching a maximum sX z ^ 1.5 — 2 and decreasing at lower redshifts |]^. Semi-analytical 
models of galaxy evolution in a hierarchically clustering cold dark matter (CDM) universe [|I| 
provide a remarkably good fit to the derived redshift dependence of the SFD. 

This scenario provides a very suitable starting point for future research. There are two 
obvious lines of investigation which urgently need to be pursued and which can be addressed 
by near-IR emission line surveys. 

1. Given the pronounced peak in the SFD between 2 ; ~ 1.5 and 2, and its remarkable 
agreement with the peak of the comoving number density of active galactic nuclei (AGNs) 
at similar redshifts and the peak in the density of damped Lya absorbers [^, this 
epoch appears to be of fundamental importance for the formation of present-day galaxies. 
However, no direct measurements of the SFD in this redshift range are available, and the 
inferred peak in the SFD follows only from the fact that both from higher and from lower 
redshifts, the SFD increases as this epoch is approached. Direct measurements of the 
SFD out to z ~ 2.5, i.e., covering the SFD peak era as well as the preceding and following 
epochs, are therefore crucially important. 

2. The SFDs at z ^ 2.75 and z ~ 4 are based on the rest frame ultraviolet (UV) emission 


of UV-bright galaxies |^. Therefore, these points are lower limits, since obscuration by 
dust will give rise to an underestimate of the derived star formation rates (SFRs) and 
since UV-obscured star forming galaxies will not be accounted for in the derived SFD. 
The colours of the HDF objects used to derive the SFD at z ^ 2.75 and 2 ; ~ 4 are 
similar to those of local UV-selected, but dusty starburst galaxies and inconsistent with 
an unreddened stellar population. Extinction corrections on the derived high- 2 ; SFRs 
of more than a factor of 10 are derived [ITl]. Obviously, the uncertainties inroduced by 
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Table 1. Redshift ranges where various bright spectral lines are available in near-IR atmo¬ 
spheric windows. 


Line 

A 

[pm] 

Z 

0.85 — 0.99 pm 

J 

1.00 — 1.34 pm 

H 

1.50 - 1.80 pm 

K 

2.03 — 2.35 pm 

Good star formation tracers 

Ha 

0.6563 

0.30-0.51 

0.52 - 1.04 

1.29 - 1.74 

2.09-2.58 

H/3 

0.4861 

0.75 - 1.04 

1.06 - 1.76 

2.09 - 2.70 

3.18-3.83 

[On] 

0.3727 

1.28 - 1.66 

1.68-2.60 

3.02-3.83 

4.45 - 5.31 

Paa 

1.8751 




0.08-0.25 

Poor star formation tracers 

Lya 

0.1215 

6.00 - 7.15 

7.23 - 10.0 

11.4-13.8 

15.7-18.3 

|0III| 

0.5007 

0.70 - 0.98 

1.00 - 1.68 

2.00-2.60 

3.05-3.69 


the presence of extinction will be strongly alleviated if a star formation tracer at longer 
wavelength can be used. 

Since Ha emission is much less strongly affected by dust than shorter wavelength tracers 
(both in continuum and in spectral lines), and since large format, high quality near-infrared 
arrays now enable emission-line surveys in the Z, J, H and R'-bands, Ha surveys will be able to 
address both of these points, and thus supply unique information on the star formation history 
of the universe out to z ~ 2.5. 


3 Use and limitations of Ho: and other tracers 


3.1 Lyo is not a good tracer of star formation 

Traditionally, surveys for high-z star forming galaxies have targeted the Lya line, which moves 
into the optical regime for z = 1.9 — 6.0. However, since Lya is resonantly scattered, even very 
small quantities of dust will effectively suppress the line @, |, ^. Indeed, despite painstaking 
efforts, Lya surveys have not revealed a population of high-z starburst galaxies [^]. Recent 
spectroscopic observations of a number of star forming z > 3 galaxies selected by the UV 
dropout technique 0 conhrm that Lya is not a good tracer of star formation in high-^ 
galaxies. These galaxies form stars at rates of ~ 10 Mq yr“^ as derived from their rest frame 
UV properties, but Lya is absent (or in absorption) in more than 50% of the cases, while in 
most of the remaining objects the line is faint [^. Thus, while Lya searches can be used to 
hnd high- 2 ; star forming galaxies ||^ 0, |^, such surveys will miss a large fraction of these 


galaxies, and cannot provide reliable SFRs for the galaxies that are detected. 


3.2 Hct as a tracer of star formation in galaxies 

The problems affecting Lya surveys are avoided by searching for Ha emission instead. The Ha 
line is not resonantly scattered and thus much less sensitive to the effects of small amounts of 
dust. In addition, the broad-band extinction at Ha (6563 A) is much less than that at Lya 
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Table 2. Conversion factors of Ha luminosity Lhq into SFR determined for various IMFs and 
empirically, and into MER, and conversion factors of various other line luminosities into Fna, 
in the sense Derived = Conversion • Observed. 


Derived 

Conversion 

Observed 

Comments 


Star formation rate ( Mq yr ^) from Lho ( L©) 

SFR 

1.27- 10-® 

Lua 

Salpeter IMF H 

SFR 

4.06 • 10-® 

Lua 

Scab IMF @ 

SFR 

1.83 ■ 10“® 

Lua 

Miller-Scalo IMF 19 

SFR 

3.41 ■ 10“® 

Lua 

empirical conversion [Q 


Metal ejection rate (MQyr ^) from Lua (L©) 

MER 

3.25 ■ 10“^° 

Lua 

not sensitive to IMF 


Lho ( L©) from luminosity in other lines ( L©) 

Lua 

2.75 


sensitive to extinction 

Lna 

9.71 

-^Paa 


Lua 

2.22 

-^[OII] 

empirical conversion 13[ , sensitive to abundance, 
excitation, IMF and extinction 


(I 215 A), RLya/^Ha = 4.28 0. Heuce an Ha survey should produce a much more reliable 
measurement of the SFD of the high-z universe. 

In Table ^ the redshift ranges where Ha and other lines are available in the near-IR windows 
are presented. Since the new 1024^ Rockwell HgCdTe “Hawaii” arrays still have at least 50% 
quantum efficiency at wavelengths as short at 0.8 /rm, the Z-band (0.85 — 0.99 /im) is included 
in this table. For Ha, the thermal background emission limits sensitive searches to redshifts 
2 ; < 2.5. Near-IR Ha surveys thus will be able to £11 in the gap between the z < 1 results of 
the CFRS and 2 ; > 2.5 results of the HDF, and cover the crucial epoch where the SFD in the 
universe reaches its maximum value. 

In interpreting Ha observations it is important to take into account that, after correcting for 
extinction at 6563 A in the rest frame, and under the usual assumption of case B recombination 
in ionization bounded, dust-free H ll regions, Ha measures only the production rate of hydrogen 
ionizing photons (A < 912A). Therefore, like the rest frame UV emission, the Ha emission 
measures only the formation rate of massive stars, and not directly the total SFR. Hence, Ha- 
derived SFRs are subject to the same uncertainties as UV-derived SFRs as far as the conversion 
from massive SFR to total SFR is concerned. This correction involves the assumption of an 


Initial Mass Function (IMF) and is both large and uncertain, as illustrated by Table pT 


However, like the UV emission, Ha emission can be used to reliably measure a metal ejection 
rate (MER), which is dominated by the same massive stars that power the Ha emission. 

A final uncertainty is introduced by the presence of dust inside the ionized regions. In 
evolved Hll regions such as are found in the disks of quiescent star forming galaxies like the 
Milky Way, absorption by dust of ionizing photons is unimportant. In (ultra)compact Hll 
regions on the other hand, dust may absorb most of the ionizing photons. In Galactic ultra¬ 
compact Hll regions, at least 50% of the ionizing radiation is absorbed by dust, and in many 
cases more than 90% [M|. In these cases hydrogen recombination line measurements, even 


at long wavelengths, will strongly underestimate the Lyman continuum flux of the ionizing 
sources. In luminous starburst galaxies, and in particular in the more extreme cases such as 
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the ultraluminous infrared galaxies (ULIRGs), star formation takes place in high density en¬ 
vironments and most of the Hll regions are compact [35|. As a resnlt, in the objects with 


the highest SFRs, Ha strongly nnderestimates the SFR. This problem cannot be solved by 
observing longer wavelength recombination lines, since the line emission is not absorbed bnt 
quenched, so that all nebnlar emission lines (and the radio free-free continnnm) will be snp- 
pressed. This phenomenon is fonnd in the prototypical ULIRG Arp 220 [^. The low Bry 
recombination line flnx of the nnclei of Arp 220 has been attribnted to a very high foregronnd 
extinction [Ay ~ 50™) However, this reasoning fails on the simple gronnds that nnclei of 
Arp 220 are directly detectable in the J, H and A'-bands so that the extinction towards 
the nnclei at these wavelengths mnst be relatively small. Fnrthermore, the low radio free-free 
emission flux density of the nuclei of Arp 220 is unaccounted for in this model. Absorption 
of ionizing photons by dust in the nuclei of Arp 220 on the other hand naturally accounts for 
all of these observations 0. 

This result has the important implication that shallow Ha surveys, probing only the most 
luminous objects are not suitable for measuring or even constraining the SFD at any redshift. 
Reliable MERs and SFDs can only be obtained if the luminosity function is probed down to 
the level where typical local Sc galaxies become detectable, i.e., at rest frame equivalent widths 
in the Ha -|- [N ll] complex EWo(Ha -|- [N ll]) > 30 A. At this rest frame equivalent width level, 
the survey depth becomes comparable to the local Ha objective prism surveys [0]. 


3.3 Use and misuse of other lines as star formation tracers 


With the exception of Lya, other hydrogen recombination lines than Ha can in principle also 
be used as star formation tracers. Since the ratios of these lines depend only weakly on physical 
consitions in the photoionized regions, only extinction plays a role in determining their relative 
strengths. Unextincted line ratios with respect to Ha for the most relevant lines are given 
in Table |3T|. The only line which can potentially produce /ory-redshift interlopers is Paa; at 
a hxed observing wavelength, all of the other lines sample higher redshifts than Ha and thus 
require much more luminous (and hence rarer) objects to yield detections. Since Paa interlopers 
always occur at low z, they are relatively easy to recognize. In contrast, optical surveys are 
strongly plagued by low- 2 ; interlopers in [O ll] and [O in] emission |^. Nevertheless, ine near-IR 
emission lines surveys it also remains necessary to spectroscopically identify the emission lines 
in all of the objects detected, since these surveys are sensitive to any emission line redshifted 
to the observing wavelength. This property has often been used to artihcially increase the 
survey volume by simultaneously considering a number of lines, assuming hxed ratios of these 
lines to Ha ll^. However, this procedure is both confusing and dangerous, since ratios of e.g.. 


[Olll] 5007 A to Ha in local galaxies cover large ranges of values |^, and depend strongly 
on excitation and abundance effects. Other than hydrogen recombination lines, only the [O ll] 
3727 A line is a useable tracer of star formation although it is considerably less well- 


behaved than the hydrogen lines, and affected by abundance and excitation effects, as well as 
extinction. More fundamentally, results based on an emission-line survey at a hxed wavelength, 
but calculated by considering simultaneously a number of diherent lines that can be detected 
at this wavelength, thus artihcially increasing the instantaneous survey volume, completely 
destroy any information on the redshift dependence of the SFD. The result is a complicated 
weighted average value over a number of redshifts, and very difficult to interpret. Thus while an 
emission-line survey at a hxed wavelength can be used to sample simultaneously volumes at a 
number of diherent redshifts in diherent lines (see Table H) , these volumes may not be combined 
to yield a “redshift-averaged” SFD, but should be treated separately, providing measurements 
at a number of diherent redshifts. 
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4 Narrow-band imaging surveys 


Blind emission-line surveys are of three types, probing different regions of parameters space 
0]: slitless spectroscopy, long-slit spectroscopy, and narrow-band imaging surveys. Because of 
the bright background in the near-IR windows, which at A < 2.2 pm is dominated by many 
narrow airglow emission lines, mainly due to OH, near-IR surveys are best conducted using 
the narrow-band imaging method. Briefly, this technique consists of deep integration in a 
narrow-band hlter (optimally chosen such that bright airglow lines are avoided) isolating a 
narrow redshift interval for a given spectral line. Objects showing a narrow-band flux excess 
when compared to a broad-band image are emission line candidates. Spectroscopic follow-up 
for securely conhrming and identifying the spectral lines and for accurately determining their 
fluxes, is required. The choice of narrow-band hlter width is determined by a trade-off between 
survey speed (which suffers if the hlters are too narrow) and equivalent width sensitivity (which 
is degraded for hlters which are too wide). 


4.1 The ESO near-IR Ha survey 


Recently, a new survey has been carried out at the ESO/MPI 2.2m telescope ^|. Based on 


the considerations of Sect. the survey was designed to emphasize depth more than area 
coverage. For the hrst time, detections were obtained. A detailed description of the survey 
parameters and observation and reduction procedures is given elsewhere 1^^. The total area 
covered is 48.0 and the total comoving volume sampled for Ha emission is 9.0 ■ 10^ Mpc^ 
for qq = 0.1 (or 3.0 • lO^h^^Mpc^ for go = 0.5), which is about 50% of the volume of the 
largest survey to date 
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1.9 • 10“^® ergs“^ cm“^ pm 
deeper. In the entire survey, two serendipitous objects were found with a narrow-band excess 
signihcant at about the da level. 

The limits implied by this survey are presented in Fig. in terms of a comoving number 
density of galaxies more luminous in Ha than a certain luminosity. The ESO survey produces 
the best currently available limits for SFRs in the range of 20 to 6 OM 0 yr“^. It is signihcant 
that the two serendipitously detected objects found in this survey have SFRs in this range. 
Locally, such SFRs are found in starburst galaxies. It is clear that in the absence of evolution 
in the Ha luminosity function no detections would be expected in any of the surveys. However, 
strong evolution is known to take place. Since both Ha and far-infrared (far-IR) emission are 
proportional to SFR, it is reasonable to assume similar evolution laws for the Ha and far-IR 
luminosity functions. The evolution of the far-IR luminosity function can be described by pure 
luminosity evolution oc (1 -|- z)^, or pure density evolution oc (1 -|- z)®’® |^. The local Ha 


. However, with an area-weighted 3cr point-source sensitivity of 
in the narrow-band frame, the ESO survey goes considerably 


luminosity function, evolved according to these evolutionary scenarios, is plotted in Fig. 
where we have continued the evolution to a redshift Zg. Inspection of Fig. ^suggests that many 
detections would be expected for Ze = 2, and a fair number for Ze = 1 as well. Since the highest 
SFRs will be underestimated by Ha (Sect. |T^) , the analysis will be limited to low luminosities. 
For Ze = 2, both luminosity and density evolution of the forms adopted here are ruled out by 
the ESO survey; in the case of density evolution, the strongest constraint is provided by the 
small-area Keck survey [Q. For Zg = 1, density evolution is fully compatible with all surveys. 


while luminosity evolution cannot be ruled out given that the Ha emission may be somewhat 
suppressed at the relevant SFRs. In summary, these results rule out evolutionary scenarios of 
the types considered here if this evolution continuous backward to Zg = 2 ; however, for Zg = 1 
the survey limits are not violated. Obviously, the same kind of agreement can be obtained by 
lowering the exponent of the evolutionary law. 
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Figure 1. Limits of the ESO near-IR Ha survey expressed as limiting comoving number 
density (vertical axis) of galaxies with Ha luminosity (top horizontal axis) or total SFR (bottom 
horizontal axis) exceeding the values indicated on the horizontal axes. SFRs were calculated 
assuming the empirical conversion factor 0 (see Table |3.1|) . The comoving number density 
limits represent 90% conhdence levels based on Poisson statistics. Ha luminosity and SFR 
limits denote 3a levels. For comparison, the results of two most important earlier surveys 
p0| , |3T| , reanalyzed in a fashion identical to the ESO survey, and considering only Ha at 
z = 2.1 — 2.4 are also shown (labeled “PD95” and “TMB96” respectively). The long-dashed 
curve represents the cumulative local Ha luminosity function p. The dotted curves show 
this cumulative luminosity function after applying pure luminosity evolution proportional to 
(1 -1- z)^ ending at a redshift Zg. Similarly, the short-dashed curves are based on the local Ha 
luminosity function assuming pure density evolution proportional to (1-1- 2 ;)® ® ending at redshift 
Zg. Evolutionary calculations for both Zg = 1 and Zg = 2 are shown. The continuous curve is 
the cumulative Ha luminosity function at ^ = 2.25 resulting from a semi-analytical calculation 
of galaxy formation and evolution in a hierarchically clustering CDM universe [M. 
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Figure 2. Implications of the z = 2.25 results of the ESO near-IR Ha survey (hlled circle) 
for the SFD (left vertical axis) and MER (right vertical axis) compared to the rest frame UV- 
derived values based on the CFRS 0 (open circles) and the HDF (open squares) and 
the local Ha survey |P] (open triangle). The vertical error bar for the ESO Ha point has been 
omitted. Star formation densities were derived from Ha luminosity densities using the empirical 
T^. Points denoted by open symbols were converted to the same scale using 
For reasons discussed in Sect. the UV-derived points 


conversion factor 
the conversion factors in Table EhT 


(open circles and squares) should be considered lower limits. 


A cumulative Ha luminosity function for = 2.25, kindly supplied by Carlton Baugh and 
based on a semi-analytical model of galaxy formation and evolution in a hierarchically clustering 
CDM universe [|^, is also shown in Fig. ^ It is remarkable that this model provides a very 
good match to the observational results. This model predicts only detections in the ESO survey, 
with SFRs of 30 — 7OM0yr“^, in complete agreement with the results. While the number of 
detections is still very small, it is clear that the ESO survey begins to approach luminosities 
and volume densities predicted by sophisticated modeling. 


4.2 Implications for star formation dut z = 2.25 


Inspection of Fig. ^ shows that only part of the Ha luminosity function at ^ = 2.25 is sampled. 
Hence, a SFD obtained by considering only the detected objects is a hrm lower limit to the 
total SFD at z = 2.25. Not unexpectedly, the resulting value of 0.014 Mq yr“^ Mpc“^ lies 
signihcantly below the value expected based on rest frame UV data . 


In order to estimate the total SFD at z = 2.25, the part of the luminosity function not repre¬ 
sented in the survey must be accounted for. The detections are quite luminous and, in terms of 
a Schechter luminosity function, likely signihcantly brighter than at z = 2.25. Consequently, 
the correction to a total SFD is large, and, since at z = 2.25 is not known, uncertain. As 
noted in Sect. |4.1| , the luminosity function for ^ = 2.25 based on the semi-analytical hierarchical 
model [Q] provides the best ht to the available Ha data, and also correctly reproduces a number 
of other observations. In the absence of a direct observational determination, this luminosity 
function is therefore the best choice for Ha emission at z = 2.25. Integration of this luminosity 











function yields a total comoving SFD at z = 2.25 of 0.12 MqYi ^ Mpc As can be seen in 
Fig.0, this value is in good agreement with the SFD expected based on the analysis by M96. 
Given the very small number of objects involved, the agreement is very encouraging. 


5 Outlook 

The above results and considerations lead to two observational recommendations. 

1. In order to provide a measurement of the SFD at a particular redshift, near-IR emission 
line surveys must probe the luminosity function down to sufficiently faint levels. Ideally, 
such surveys should aim to determine the luminosity function of star forming galaxies at 
the relevant redshifts directly, removing the need of assuming a luminosity function. The 
new generation of 8 m class telescopes equipped with large-format near-IR array cameras 
will provide the required sensitivity and survey power. By carrying out narrow-band 
searches at a number of discrete redshift intervals, selected to avoid bright airglow lines 
in the narrow-band hlters, near-IR Ha surveys will trace the cosmic star formation history 
out to z ^ 2.5, including the important = 1.5 — 2 era where the cosmic star formation 
density is believed to peak. Higher redshifts will be available by using H/3 or [O ll]. These 
surveys will be much less sensitive to uncertainties and biases introduced by extinction 
than rest frame UV based searches. 

2. The brightest part of the galaxy luminosity function cannot be probed using nebular emis¬ 
sion line surveys. Here total SFRs must be derived from the rest frame far-IR continuum. 
At low z, the brightest, optically obscured part of the luminosity function accounts for 
only a small fraction of the total SFD. However, if strong luminosity evolution takes 
place, this fraction may become signihcant at higher redshift. Blind submillimetre sur¬ 
veys (sampling the far-IR continuum of high- 2 ; galaxies) are required to address this point 
observationally. The FIRST mission will play a major role in this held. In the immediate 
future, a hrst assessment of the importance of dusty, very luminous star forming galax¬ 
ies at high 2 ; will be made with the SCUBA submillimetre camera at the James Clerk 
Maxwell Telescope. 
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